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We present finite bias measurements on a quantum dot coupled capacitively to a quantum point 
contact used as a charge detector. The transconductance signal measured in the quantum point 
contact at finite dot bias shows structure which allows us to determine the time-averaged charge on 
the dot in the non-blockaded regime and to estimate the coupling of the dot to the leads. 

PACS numbers: 73.21.La, 73.23.Hk 



I. INTRODUCTION 

The charge state of a quantum dot can be read out 
using a nearby quantum point contact (QPC) as a 
detectori Real-time readout has been recently demon- 
strated using radiofrcquency single electron transistors^ 
or QPCs^* Together with a spin-charge conversion 
mechanism, this is being considered a candidate for a 
qubit readout scheme^ in a future quantum computing 
device based on coupled quantum dots£ 

A valuable piece of information about a quantum dot's 
characteristics is the knowledge about its coupling to 
each of the reservoirs. This coupling is determined by 
the electrostatic barrier forming the constriction and the 
wave function overlap leading to tunneling. The latter 
may strongly depend on the quantum state under con- 
sideration in the dot, which means that the quantum me- 
chanical tunnel coupling has to be determined for each 
state individually. In the Coulomb blockade regime, for 
the case of single- level transport, a fit of a transport 
peak in the Coulomb blockade regime enables one to ex- 
tract an effective tunnel coupling r o ff = TsFd/^s + Td) 
where Ts and Td are the couplings to the source and the 
drain reservoir, respectively. However, the ratio Ts/Ed 
remains unknown. 

In earlier work on finite bias transport^, e.g. in car- 
bon nanotubes^ it was shown that in principle, the in- 
dividual coupling to both leads of the dot can be ex- 
tracted from the current amplitudes at positive and neg- 
ative bias. Since these considerations make use of a spin 
blockade effect, however, they apply only to the case of 
spin degenerate states in the dot, and require an even 
number of electrons on the dot. In addition, they rely 
on the absence of cotunncling. Recently, a method was 
presented^ to measure the coupling of a dot to its reser- 
voirs for a three-terminal quantum dot. 

In the following, we present finite bias measurements 
of transport through a quantum dot, complemented 
by simultaneous measurements of the conductance of a 
nearby, electrostatically coupled QPC used as a charge 
detector. While the latter allows us to determine the 
time-averaged charge on the quantum dot even in the 



non-blockaded regime, the combination of both methods 
makes it possible to provide qualitative estimates for the 
quantum mechanical coupling of the dot's energy levels 
to each of the two reservoirs. 



II. EXPERIMENTAL SETUP 

The sample (see Fig. QIa)) was fabricated using surface 
probe lithograph} ! 9 : 10 : 11 on a GaAs/Al .3Gao.7As het- 
erostructure, containing a two-dimensional electron gas 
(2DEG) 34 nm below the surface as well as a backgate 
(BG) 1400 nm below the 2DEG. The unstructured 2DEG 
had a mobility of (3.5 ±0.5) ■ 10 5 cm 2 /Vs and an electron 
density 4.6 ■ 10 11 cm -2 at a BG voltage Vbg = -0.5 V at 
T = 4.2K. 

All measurements were performed in a dilution refrig- 
erator at a base temperature of 80 mK. 

Negative voltages were applied to the surrounding- 
gates (Gl, G2, Sqpc, Dqpc, the latter two also con- 
taining the charge detection circuit; see Fig. |HV)). and 
to the BG, to reduce the charge on the dot and close its 
tunnel barriers. A voltage applied to gate P was used to 
tune the detector QPC to a regime where it is sensitive 
to the charge on the dot. The QD bias voltage was ap- 
plied symmetrically (with respect to ground) across the 
dot between source (Sqd) and drain (Dqd). 

Due to the electrostatic coupling of the QPC to the 
dot, a change in the dot's charge leads to a modifi- 
cation in the QPC's confining potential, resulting in a 
change of its conductance^ The latter was measured 
by applying a dc voltage and measuring the resulting 
current. Each additional electron on the dot leads via 
electrostatic interaction to a shift of the QPC conduc- 
tance's dependence on gate voltage: / qP c,N(V^atc) = 
^qpc,N+i(Vj»atc + AV), where N is the number of elec- 
trons on the dot. If we assume that the QPC is tuned 
to a regime between two conductance plateaux with an 
approximately constant derivative dlqpc/ dVg&tc, then we 
can subtract Iqp C ,7v=const(^gato) as a background and get 
a signal proportional to the additional charge on the dot: 
Q oc Iqpc(Vgate) - -fq P c,JV=const ( Vg a to ) ■ This technique 
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FIG. 1: (a) AFM micrograph of the structure with designa- 
tions of gates: source (S) and drain (D) of the quantum dot 
(QD) and the quantum point contact (QPC) used as a charge 
detector; lateral gates Gl and G2 to control the coupling of 
the dots to the reservoirs; plunger gate (P) to tune the QPC 
detector. Only the part of the circuit related to the readout 
functionality is included. Vq pc consists of a dc and a small 
ac component (Vq pc ,ac < 5/xV). (b) Example measurement 
of the current through the dot. For this low bias measure- 
ments, the voltage across the dot was Vbias, dot = 10 yN . (c) 
Simultaneous measurement of the conductance through the 
QPC, where each step corresponds to a change of the dot's 
charge by one electron, (d) Simultaneous measurement of the 
transconductance dl qpc / dV^pc. 



has recently been applied to investigate the charging be- 
haviour of a double quantum dot^ 

In our measurements, in order to avoid the influence of 
nonmonotonuous drift in I qpc , we measured the transcon- 
ductance dl qpc / dV q p C as an additional quantity, using a 
lock-in setup at a frequency of / = 31 Hz. The volt- 
age V^p C was composed of a constant dc voltage, used to 
tune both the dot's and the QPC's chemical potential, 
and a small ac voltage used to periodically change the 
chemical potential inside the dot by a small amount (of 
the order of k-oT). The bias Vbi a s,qpc across the QPC 
was applied symmetrically with respect to Vq pc , in or- 
der to minimize its influence on the chemical potential 
inside the dot. Figure ^a) illustrates the QPC related 
part of the circuit diagram. Figures ^b)-(d) show the 
correlations between the three quantities in a simultane- 
ous measurement: at the position of a Coulomb block- 
ade peak in the current through the dot, a kink appears 
in the QPC's conductance, corresponding to a change 
in the dot's charge by one. In the transconductance 
dI qpc /dV qpc , a dip is observed. dI qpc /dV qpc , being a 
derivative of the QPC current, differs from dl qpc / dVQ2 
only by a constant background and a factor given by the 
ratio of the corresponding lever arms, o;G2/aqpc- This ra- 
tio should depend only weakly on gate voltages. One can 
therefore obtain the time-averaged charge on the dot by 
integrating the measured value dI qpc /dV qpc with respect 
to Vq2, (after subtracting a constant and a weak linear 
background to compensate for the gate voltage depen- 
dence of ciG2/ciqpc) and normalizing the resulting steps 
to unity. 
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FIG. 2: (Color online) (a) Finite bias transport measurement 
through the dot at a magnetic field of B = 0.1 T. Dark re- 
gions represent larger positive [blue] or negative [red] differ- 
ential conductance, (b) Corresponding measurement of the 
transconductance. Dark regions represent larger values, (c) 
Single transconductance trace for low bias (Vbias = 20 //V, 
see upper horizontal line in (b)). (d) Single transconductance 
trace for Vbias = — 0.2mV (see dashed horizontal line in (b)). 
(e) Averaged transconductance. Averaging was performed 
over a set of bias voltages — 1.6mV < Vbias < — 0.6mV. Indi- 
vidual traces were laterally shifted with respect to each other 
so that the two lines visible in the lower left part of (b) (par- 
allel to the diamond edges in (a)) yield two distinct sharp 
peaks. The region delimited by dotted lines in (b) marks the 
range over which averaging took place. Furthermore, a con- 
stant and a linear background were subtracted, (f) Averaged 
transconductance integrated with respect to Vg2- 

III. TIME AVERAGED CHARGE DETECTION 

In Fig. |2t a), we present finite bias measurements of the 
dot's conductance. From these, estimates for the charg- 
ing energy E c w 1.6 me V and the mean level spacing 
A w 300 /ieV were extracted. Fig. |2Ib) shows a simulta- 
neous measurement of the QPC's transconductance. 

In the transconductance plot, diagonal lines with the 
same slope as the diamond boundaries in the dot conduc- 
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tance plot are observed, marking a change of the dot's 
time-averaged charge. The lines correspond to the align- 
ment of an energy level in the dot with either the source 
(negative slope) or the drain (positive slope) reservoir, 
while their intensity contains information about the mag- 
nitude of the change in charge. 

Figure |5{c) shows the transconductance at low bias. 
Dips in the signal at the gate voltages of the Coulomb 
blockade peaks correspond to the change in electronic 
charge by one elementary charge each. For a finite bias 
Vbms = — 0.2mV (Fig. Old)), the peaks split, illustrating 
that the time-averaged charge on the dot changes in steps 
smaller than the unit charge e. 

These steps in charge can be directly visualized by in- 
tegrating over the transconductance (Fig. Ele)-(f)). In 
order to improve the precision of the charge determina- 
tion, the transconductance was averaged over a finite bias 
range prior to numerical integration. A constant (cor- 
responding to the direct coupling between the in-planc 
gates used for ac excitation and the QPC) and a linear 
term (corresponding to the gate voltage dependence of 
c*G2/a!qpc) were then subtracted. To conserve the steep- 
ness of the steps near Vg2 ~ — 240 mV, the individual 
traces were laterally shifted with respect to each other 
before averaging. As a consequence, the step correspond- 
ing to the second Coulomb peak is broadened, since for 
this peak the slope of the corresponding lines in Fig. |5f b) 
is different. 

The large steps marked by dashed horizontal lines in 
Fig- Elf); corresponding to a change in charge by one 
elementary charge each, have almost identical height, as 
one would expect. Charge rearrangements during the 
measurements, one of which is visible at Vg2 ~ 235 mV 
in Fig. 2(a), might contribute to errors in the numerical 
integration and thus may lead to small differences in total 
step height. 

To verify the usefulness of the integrated signal as a 
measure for the mean charge on the dot, we compared 
the step size related to a single electronic charge for sev- 
eral Coulomb peaks (see Fig. |2Jf)) and for different bias 
ranges between zero and the height of a Coulomb dia- 
mond. The difference in step height did not exceed 10 
percent if a sufficient number of traces was taken into 
account to reduce the statistical error. Even in the high- 
bias regime 1.2 mV < Vbias < 1.6 mV near the border 
of the Coulomb-blockaded bias range, where the lines in 
the transconductance seem to gradually disappear, the 
integrated transconductance signal gives a good measure 
of the mean charge, which means that line broadening 
compensates for the lower amplitude. 



IV. ESTIMATION OF DOT-LEAD COUPLING: 
NON-DEGENERATE CASE 

From the information about the mean charge in the 
non-blockaded region, it is possible to extract the cou- 
pling of individual states to each lead separately by us- 



ing a rate equation approach*^ if a single energy level lies 
within the (thermally broadened) bias range between the 
Fermi level of both leads, the mean dwell time of the Ath 
electron in the dot is determined by the relative value of 
the couplings rs,jv and Td.at to both leads. The Flat's 
account for the height of the tunnel barriers i, the wave 
function overlap between the dot and lead, and the den- 
sity of states inside the leads, which we assume to be 
constant. The occupation probability of the energy level 
A then becomes 



P(N) 



rs,w/s(/- t A r ) + To, at jb (mat) 



N 



(1) 



where / stands for the energy distribution in the leads 
(which we assume to be the Fermi distribution) and the 
indices S and D stand for source and drain, respectively. 
For a bias voltage | Vbias | = \Ep,s — #f,d| ^> k^T and a 
chemical potential /in not inside the thermally broadened 
ranges around the Fermi energies of both leads, we have 
fs(HN) = I? /dO^aO = 0; an d this expression becomes 
P+(A) = r SiA r/(rs,A/ + r DiA r) for E F . S > £ FiD (positive 
bias) and P~(N) = r DiA ,/(r s ,A/ + r DiA r) for E F . S < £ f ,d 
(negative bias). 

Since the non-integer (fluctuating) part of the mean 
charge on the dot in the level A is given by eP ± , this 
connects the integrated transconductance value to the 
ratio T s . n /T u ,n- 



T S , N P+(N) 



P H 



(2) 



r D jv p-(n) i-p+- 

Together with the expression for the transport current 



rs,7vrD,jv 



(3) 



it is possible to determine the numerical values of Ts } n 
and Fd,at. 

The analysis of the transconductance signal becomes 
more involved if more energy levels inside the quantum 
dot are contributing to transport. It is thus preferable to 
start the analysis in the regime of single- level transport, 
i.e. Vbias < (A — ksT). In addition, the maxima in 
the transconductance signal associated with source and 
drain alignment of the level must be clearly separable, 
i.e. Vbias > k^T. Those two requirements limit the bias 
range over which the integration of the transconductance 
peak can be performed on a given data set, especially 
when excited states are present near the ground state. 
To apply the method described above, we have to verify 
that there exists a finite bias range where these conditions 
are met. 

We analyzed the low finite bias regime around three 
conductance peaks in a range where the system was most 
stable, for five different magnetic fields applied, to change 
the shape of the dot's wave functions and their overlap 
with the leads. In our analysis, we encounter two quali- 
tatively different situations in the low bias transconduc- 
tance regime: No branching at zero bias (this section): 
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single line with constant slope crossing the zero bias line, 
possibly branching at finite bias. Branching at zero bias 
(degenerate case, discussed in the next section). 

An example of the first case is encountered, e.g., at the 
second peak in Fig. Et a )/(b) (sec detail in Fig. EI a )), 
in the low bias regime: near Vbi as = 0, we observe a 
broadened line with negative slope in the transconduc- 
tance, while a similarly broadened, but comparatively 
weak maximum is measured in the transport current 
(Fig. Ha)). 

Figure Elb) shows a schematic representation of Fig. 
OUa), illustrating which line corresponds to the alignment 
of the dot's chemical potential with source (S) and drain 
(D). Figures EJc)-(f) illustrate how coupling and level 
alignment determine the mean electron number on the 
dot. Note that in Fig. Etc), an asymmetrically coupled 
level can become nearly 100 percent occupied even in the 
non-blockaded regime. 

The observation near Vbias = is compatible with the 
presence of a single energy level strongly coupled to the 
source reservoir (leading to enhanced broadening, com- 
pared to a Coulomb blockade peak width of purely ther- 
mal origin), but extremely weakly coupled to the drain 
contact. This would explain the observed low zero-bias 
transport current and the fact that within our measure- 
ment resolution, no line with positive slope correspond- 
ing to an alignment to the drain contact is observed in 
the transconductance. At finite negative bias, i.e. once 
a more symmetrically coupled excited state becomes ac- 
cessible (see finite bias region around middle peak in Fig. 
EJa)), transport is strongly enhanced. 

Solving equations ^ and El f° r Ts.jv and Td.at yields 
an upper bound for the ratio of the coupling constants 
^d.n /Ts.n < 1/50. From the measured value of the 
transport current in the finite bias regime (averaged over 
positive and negative bias), we obtain an estimate for 
Ts,ivrD,A»y (rs,jv + Td.at), so that both are uniquely de- 
termined within the errors mainly resulting from charge 
measurement uncertainties. In fact, for this very asym- 
metric case, we might assume that level broadening is 
caused entirely by the more strongly coupled source lead. 
The full width at half maximum of the Coulomb peak 
is of the order iufwhm = 130 [ieV < A, yielding a 
value for Tg^ which suggest an even stronger asymmetry 

r D ,iv/rs,jv < 1/1000. 

At finite bias, a splitting of the line in the transcon- 
ductance is observed, suggesting that an excited state of 
the dot more symmetrically coupled to both leads gov- 
erns its mean occupation. Inserting this into a numerical 
simulation (see below) also yields qualitative agreement 
with the observed finite bias dot current, which shows a 
larger step change at the positive than at the negative 
bias corresponding to the excitation energy. 
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FIG. 3: (a) Detail from Fig. ffib). showing a broadened 
line with negative slope at zero bias, and a splitting at finite 
bias, (b) Illustration of the connection between the coupling 
of a single state (G) to source (S) and drain (D) and the 
mean charge on the dot. This representation was generated 
using a simulation based on a rate equation approach,— ^ and 
is compatible with our experimental observations (with the 
exception of line broadening due to strong coupling, which 
is not included in the theory). Our model for this case in- 
cludes two single-particle levels: an asymmetrically coupled 
state accounting for the imbalance of line weight observed 
at positive and low negative bias (where the line with posi- 
tive slope, corresponding to alignment of G with D, is practi- 
cally absent), and a higher energy level H with more symmet- 
ric coupling, responsible for the branching at finite negative 
bias. We used coupling asymmetries Fd,g /Ts.g = 1/50 for 
the level G, and TD.H/rs.H = 1/4 for the level H. In addi- 
ton, we set Td.g + Ts,g = Td,h + Fs,h- The diagonal line 
(G-S) with negative slope corresponds to the alignment of 
a single energy level (G) with the source (S) contact, while 
the line (H-D) with positive slope in the negative bias range 
marks an alignment of an excited state (H) with the drain (D) 
contact. Letters in brackets (c)-(f) located in the blockaded 
(left and right) and non-blockaded (top and bottom) regions 
refer to the corresponding figures. These illustrate how the 
time-averaged electron number depends on bias and chemi- 
cal potential of the dot for the case of an asymmetric (here: 
P + > |) coupling. 



V. ESTIMATION OF DOT-LEAD COUPLING: 
DEGENERATE CASE 

In a number of cases the lines in the transconductance 
plot branch at zero bias (leftmost peak in Fig. Elb), see 
detail in Fig. Ufa)). This splitting is not compatible with 
the participation of only one individual single-particle 
level in transport. At finite positive bias, a further split- 
ting occurs and a line with opposite bias dominates, sug- 
gesting predominant coupling towards the source contact. 

From our transconductance measurements, the mean 
occupation was determined by integrating over a volt- 
age range in Vq2 in the transconductance plot, summing 
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FIG. 4: (a) Detail from Fig. ISTb). showing discontinu- 
ity/splitting of lines in the transconductance signal at zero 
bias, (b) Schematic representation of (a) generated by a nu- 
merical simulations. To explain the observed line splitting, a 
simplified model spectrum with four states is used, consist- 
ing of two pairs of (quasi-)degenerate levels. Small letters 
in brackets refer to Figures (c)-(e), where the spectrum is de- 
picted for different values of gate and bias voltage, illustrating 
how the interplay of several levels with different couplings to 
the leads influences the time-averaged occupation of the dot. 
The labeled lines in (b) refer to the alignment of a level (L, G, 
or H/H', see Fig. (c)) with either source (S) or drain (D). (c) 
For low positive bias, charge passes through the degenerate 
levels G and L. Due to the asymmetric coupling, the mean 
charge remains low until alignment with the drain contact, 
(d) For small negative bias, the asymmetry in the coupling 
of levels L and G leads to a finite mean charge on the dot 
as soon as they align with the drain's chemical potential (see 
line G-D in Fig. jl^b)). (e) In this parameter range, the mean 
charge reaches a high level, due to the asymmetric coupling 
of the levels H and H' that trap charge entering through the 
source contact. 



over a range of bias voltages to reduce the statistical er- 
rors. Normalization was performed by integrating over 
the whole non-blockaded gate voltage range, summing 
over the same bias range. 

To elucidate the origin of the low bias results (depicted 
schematically in Fig. Efb)), we performed analytical cal- 
culations using a rate equation approach based on the 
framework of Bccnakker's theory of sequential electron 
tunneling. 13 

We assumed that a model involving at least two quasi- 
degenerate states might describe our findings, and there- 
fore extended eauations lll3l to the case of two degenerate 



levels, involving four coupling values r ?i 



S,D for 



source and drain; j ~ 1, 2). 

To extract a meaningful result despite the large relative 
errors in the charge input values, we express the l\j in 

terms of a global multiplicative factor r = ^/Ili , 



< r S:i 
r D ., 



rs.iTD, 



the symmetries of each of the two states Si 

and the relative weight of the two states V = 

expressed as the ratio of the geometrical averages of their 
couplings. The resulting simple expressions for the I\ j 
become: r s ,i - T VSx, r D ,i = T V/Si, T s . 2 = T a S 2 /V, 

r D , 2 - r /(V5 2 ). 

Using the four charge and transport values at low 
positive and negative bias as input values, the expres- 
sions for current and occupation can be solved for the 
four new parameters. The numerical values obtained are 
fiT = (2.0±0.7)/icV, V = 1.23 ±0.50, Si = 1.10 ±0.15, 
5*2 = 0.22 ±0.12. This means that one of the two lev- 
els has approximately symmetric coupling, while the sec- 
ond is predominantly coupled to the drain contact. This 
state must therefore have a very asymmetric distribution 
of the wave function amplitudes to the two reservoir con- 
nections, since we learn from the slopes observed in Fig. 
I2b) that the dot as a whole predominantly couples to the 
source reservoir, suggesting an asymmetry in the tunnel 
barriers themselves. Since the numbers show that the 
two states have different symmetry properties, the most 
obvious possibility of a spin pair has to be excluded in 
this case. 

The method described above can be extended to deter- 
mine the coupling properties of higher energy states by 
incrementally solving the more intricate rate equations 
involving a larger number of levels, using measurements 
at higher bias and the previously calculated r,-^ values as 
an input. The major advantage of this procedure is that 
it correlates data measured at a constant value of the 
controlling gate, changing only the bias of the dot: due 
to spectral scrambling, the information obtained in this 
way is rarely accessible by comparing successive Coulomb 
peaks. 

Wc pcrfomcd this analysis using numerical calcula- 
tions based on the same aforementioned rate equation 
approach. Even though the input values include already 
a certain error, some qualitative statements can be de- 
rived: the most important one concerns the feature vis- 
ible in Fig. Ufa), showing a change in slope in the right 
positive bias line occuring at finite bias. The fact that the 
line with positive slope is discontinued is not compatible 
with the presence of a single excited state. Again, two 
quasi-degenerate states are necessary to explain the sce- 
nario observed. The coupling of these two states has to 
be asymmetric with predominant coupling to the source 
contact. The measurements are compatible with the two 
states having identical coupling values, leaving open the 
possibility of a spin pair. Fig. EJb) shows a numerical 
simulation using the values T^.j from the analytical cal- 
culations above for the two lower states (responsible for 
the low bias behaviour). The remaining Ti_j were deter- 
mined numerically by iteratively comparing experiment 
and simulation. 

The resolution of our charge detector is limited: all the 
transconductance measurements presented above were 
performed in a regime in which the QPC's conductance 
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remained between 2e 2 /h and Ae 2 /h. Slight changes in the 
lithographic pattern for future structures will permit us 
to approach the tunneling regime, resulting in a tenfold 
increase in charge sensitivity and a more precise deter- 
mination of the mean charge and the coupling strengths. 

VI. CONCLUSIONS 

We have demonstrated how a QPC charge detector 
for a nearby quantum dot can be used to extract infor- 



mation about the average charge in the non-blockaded 
finite-bias regime. Together with the conductance mea- 
surement through the dot itself, this allows us to extract 
qualitative and quantitative information about the cou- 
pling of single energy levels to the source and drain reser- 
voirs. 

Financial support from the Schweizerischer National- 
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